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We investigate the modification induced on the nucleon-nucleon cross section by 
screening of the interaction inside nuclear matter. The interaction is described by the 
relativistic one boson exchange of a, uj and 7r mesons. Medium effects enter as loop cor- 
rections to the meson propagators and are characterized not only by density, but also by 
temperature and velocity with respect to the center of mass of the collision. The cross 
section displays peaks for some specific values of the velocity, corresponding to the excita- 
tion of zero-sound modes of the longitudinal u> field. The enhancement factor amounts to 
about 1.5 under reasonable conditions. It increases with density and is reduced at finite 
temperature. These findings may have verifiable consequences on the angular dependence 
^ ' of the measurable secondary particle distributions. 



1 Introduction 

Dissipative effects play an important role in explaining the shapes of the angular distributions 
and energy spectra of the secondary particles produced in a relativistic heavy ion collision. For 
this reason, the cross section is an important ingredient in all kinds of numerical simulations 
of these collisions, be it in the molecular dynamics approach 0]j kinetic theory approaches 
[§, or the hydrodynamical models where it enters through the transport coefficients ||. 
It has become obvious in the last years that using the free value of the cross section does not 
well reproduce the available experimental data; rather, the medium may affect the nucleon- 
nucleon cross section by modifying the strength of the interaction and the available phase 
space. 

Several attempts to calculate medium corrections to the cross section already exist. The 
first point to precise is in what type of background matter does the collision take place. Some 
authors like Faessler et al. || or Fuchs et al. ^ consider that the background consists of a 
momentum distribution in two well separated Fermi spheres. This corresponds to the initial 
stage of a heavy ion collision and is far out of equilibrium. Other authors like Machleidt et 
al. H or Aim et al. [[| consider that the background is nuclear matter, supposedly of infinite 
extension, and not far from equilibrium. This rather corresponds to the intermediate stage 
of the collision after a fireball was formed. This latter case is what we are interested in here. 
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Medium effects enter at several levels of the calculation of the cross section. Older eval- 
uations only took into account the shifting of the nuclear mass to an effective smaller value 
inside nuclear matter. Then, in the last few years were performed a number of calculations 
with the Brueckner formalism £j pi]] . In this approach, an infinite series of ladder diagrams 
are resummed. The nucleon propagators, i.e. the uprights of the ladder, are dressed. The 
medium effects enter as Pauli blocking factors on the intermediate states in the inner nucleon 
lines. The meson lines are the free ones, corresponding to an unscreened interaction. 

On the other hand, we would like to discuss in this letter the influence of the screening of 
the interaction by the surrounding nuclear matter. In a model where the N-N interaction is 
described by relativistic meson exchanges, we take into account the dressing of the meson lines 
in the first Born approximation. A detailed study of the density and temperature dependence 
of the total and differential elastic p-p cross section was performed in our first paper [11]. It 
was assumed there that the background nuclear matter was at rest with respect to the center 
of mass of the collision. The concern of our present work is to release this assumption. We 
refer the reader to [11] for details on the formalism and only recall here the main formulaes. 



2 The model 

We are working with the standard quantum hadrodynamics lagrangian with a, u and ir meson 
exchange. This is the minimal model which permits to reproduce the available experimental 
data on p-p scattering to a satisfactory level of accuracy. The interaction lagrangian reads 

C = g a ip<Tip + 9u> Ip^LOpIp + g^lp^TTTlp (1) 

The differential elastic cross section was calculated in the Born approximation. In the center 
of mass of the collision, 

da/dSl = l/(64vr 2 s)|A4| 2 (2) 
where s is the Mandelstam variable s = 4E 2 m and the transition matrix is given by 

M= J2 (P3T m Ui)G$t(U 4 r n U2) - exchange (3) 

m,n=(T,a;,7r 

In this expression the f/'s are the spinors defining the incoming [U\, U2) and outgoing (C/3, U4) 
nucleon states. The indices i = 1, 2, 3, 4 carried by the spinors stand for the momenta pf and 
the spin-isospin quantum numbers Sj, U. Crfi" is the propagator matrix which depends on 
the momentum transfer k^ = p^ — p% . In free space, it reduces to the standard form 

Gfi = Diag f ~% U ~T ] > V^-^ ) W 

\ k z — m z k z — m z fe z — m z j 

r is the coupling matrix T = (—ig^^, ig a , 9^75'?) • in the actual calculation of the transition 
matrix, all coupling constants are multiplied by form factors (A^ — m 2 a )/(A 2 n — k 2 ). The 
exchange term is obtained by interverting indices 3 and 4. In the following we consider 
only the spin-averaged cross section |A4| 2 = (1/4) £) a . \A4\ 2 . a vv and a nn are equal in our 
approximation where the Coulomb interaction was substracted and the neutron and proton 
masses are equal. We do not show a np since our simple model does not include p meson 
exchange. 

In paper I we performed a fit of the free elastic p-p cross section using the available 
experimental data [jl2|l . We kept constant the masses of the mesons and the pion-nucleon 
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coupling (^=550 MeV, m w =783 MeV, m 7r =550 MeV, #£/47r=14.4) and fitted the coupling 
constants g a , g^ and the cutoffs A CT , A w , A,,-. Several sets of parameters were found to 
reproduce the data to a satisfactory level of accuracy. It was checked that the in-medium 
behavior of the cross section was the same for all choices of the parameter set. The calculations 
are performed here with parameter set A of paper I: g a = 3.80, g u = 9.31, A CT = 1298.8, 
A w = 1240.5, A,, = 362.1. 

Inside the medium, the momenta and masses of the incoming and outgoing nucleons must 
be replaced by their effective counterparts. In the Hartree approximation, — » (pf )* = pf — 
g^ <u^>, m — > m* = m — g a <a>. The propagator matrix must also be modified in order 
to take into account that the N-N interaction is mediated by dressed mesons. The dispersion 
relations and propagators of the mesons at finite temperature were obtained from a linear 
kinetic analysis of perturbations around the Hartree ground state jy|]. The results of this 
method coincide at T=0 with those of the one- loop approximation |l4| . The calculation takes 
into account renormalized vacuum polarization contributions which are crucial in obtaining 
a physically reasonable behavior of the propagation modes 13]. The propagator matrix so 
obtained has the form ][ll|, 14, 15 1. 
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In these expressions we defined the quadrivector rf and the projectors onto the longitudinal 
(A^) and transverse {T^ u ) modes as functions of the 4-velocity of the fluid as follows: 



k.u 



rf rf 
rj 2 



n /JtV 



Wk h 



y-jf-ly-jlS 

rf 2 



(7) 



The polarizations H(k) of the meson fields are functions of the thermodynamical state (density 
and temperature) of the matter and on the transferred 4-momentum k^ = (to, q). There is a 
mixing between the a and longitudinal-u; propagation modes. The 7r propagation is decoupled 
from the dynamics of the other mesons in symmetric nuclear matter. We use the notation 



IT 



n„ 



-A TT^ 



u. 



-Vnr/2 



(8) 



for the mixing, longitudinal and transverse parts. The explicit expression for the polarizations 
and a study of the dispersion relations can found in jll|, 13]. Besides the a and longitudinal 
and transverse u branches, there is a zero sound branch in the mixed <j-oj mode. 

The denominator of the a, longitudinal part of the u and mixing propagators is given by 
the mixed g-oj dispersion relation 



Dv-uLik) = (k 2 -m 2 a + glU a )(k 2 - ml + glU^) + g a gjl 



(9) 
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The polarizations are complex. In the T = limit, the imaginary part of the polarizations 
vanish outside the region defined by 



V \pf ~ q) 2 + ml - \J 'pp + ml < u < y/ (p F + q) 2 + ml - \J p 2 F + ml (10) 

(where pf is the Fermi momentum and m* the effective mass of the nucleon inside the 
medium). Inside this region, the propagation modes are damped by the decay into particle- 
hole pairs. At finite temperature the imaginary parts are finite for any space-like mode. 

The dependence on the velocity of the medium comes through the quadrivector rf". In 
the referential (Rl) where the fluid is at rest, = (1,0) and the momentum transfer is 
k^ = (oj, q) 

In the center of mass of the collision (referential (R2)), the background fluid moves with 
a velocity u M = (7, 72?) parametrized by 

v = (v sin a sin (p, v sin a cos cp, v cos a) (11) 

The momentum transfers for the direct (&31) and exchange (fen) diagrams are 

= (0,0,psm9,p(l - cos 9)) ; fe£ = (0, 0, -psin0,p(l + cos<9)) (12) 

with p = \J E'i m + m'l and 9 is the cm. scattering angle. 

One goes from (Rl) to (R2) by performing a boost of four-velocity (7, jv) and a rotation. 
In the referential (Rl), 

(ujQ,q e ) ; kl^Rl) = (w e) ?e) 

^vp[± sin a cos ip sin 9 — cos a(l =F cos 9)] (13) 
[wf +2p 2 (lTcos#)] 1 / 2 



k^(Rl) = 

q o = 



3 In medium cross section for finite relative velocity 

When the rest frame of the background nuclear matter is moving with respect to the center of 
mass of the colliding particles, the relative velocity will enter the expression of the cross section 
through the meson propagators. There is an explicit dependence coming from the projectors 
, T^ u (see Eq. (^)). The projectors act by attributing velocity-dependent weights to the 
various components of the currents J like e.g. J£ v = (Uzj^Ui) {U^Uz). There is also 
an indirect dependence in the velocity in the expression of the polarizations IL(u) ( e^,q ( e^) 

coming from the boost transformation Eq. (13). In contrast to the zero velocity case where 
the momentum transfer entering the expression of the polarizations was confined to the 
axis u=0, k^ = {uj ( e\,q,e\) now spans all the spacelike region. We will consider three cases: 

(i) v = (0,0,0 

For moderate values of the velocity along the axis defined by the incident particles (v = v z ), 
the cross section is slightly enhanced at low densities and reduced at high densities. When 
some temperature is introduced, the low density enhancement disappears and the high density 
reduction persists. At highly relativistic v z velocities ( r y z ~ 5), the cross section displays a 
series of peaks, of which three are clearly visible. The height of the peaks increases with 
density. The effect of finite temperature is to gradually wash out these peaks. The peaks are 
superposed on an additional increase of the cross section, which stabilizes to a large constant 
value for large j z . This additional enhancement persists when introducing finite temperature. 
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(a) v = (o,v y ,o) 

For low values of the velocity inside the reaction plane and orthogonal to the beam axis 
(v = v y ), the cross section is first approximately constant for densities up to n = uq and 
somewhat reduced at high density when increasing v y . Further increasing v y , one meets a peak 
structure (see Fig. 1) superposed on a smooth increase similar to the one already observed 
in the case v = v z , but occurring for more physical values of the velocity. There is first a 
sharp peak for a characterictic value v y immediately followed by one or two smoother peaks. 
Further increasing the velocity, the cross section slowly decreases and stabilizes to a constant 
value (Jmed/ & free( v y —* c) somewhat larger than (T me( i/crf ree (v = 0) and approximately equal 
to 1. We were able to check that the peak structure originates from the velocity dependence 
of the polarizations while the residual temperature-independent enhancement at v y — ► c is 
due to the velocity dependence of the projectors A^, T^ v . The same holds for the case v//z. 

(Hi) v = (^,0,0) 

No appreciable modification is observed for velocities orthogonal to the reaction plane (v = 
v x ) up to v x ~ 0.8c. On the other hand, the total cross section starts to increase above this 
value and continues to grow indefinitely as ^ x — > oo. 

The angular dependence of the differential cross section deviates appreciably from its 
behavior when the fluid is at rest. For each value of the velocity, there is a favoured angle, 
but this effect is most appreciable for small angles, so that finite velocities also contribute to 
make the cross section more forward/backward peaked with respect to its vacuum value. 

We concentrate on the behaviour with v y which is the case most liable to lead to observ- 
able consequences since the peaks appear for reasonable values of the relative velocity. The 
characteristic velocity v y at which the sharp peak takes place is density dependent, e.g. one 
has v y ~ 0.41c for n=riQ and v y ~ 0.57c for u=2uq (see Fig. la). The height of the peak 
increases with increasing density and is reduced by increasing the temperature (see Fig. lb). 
At T=10 MeV, a broad bump is still clearly visible around v y , at high temperature T ~ 50 
- 100 MeV it has completely disappeared. 

The peak structure is to be related to the excitation of the zero sound modes. These 
modes arise from the mixed a - longitudinal uj part of the dispersion relation; they were 
studied in previous work [13, 15]. A contour plot of the mixed dispersion relation ^-^(w, q) 



in the q — uj plane of transferred momentum is shown in Fig. 2 for T=0. h=uq. The level 
D a -ujL(u,q) = gives the location of the zero-sound branch. For not too high values of 
the energy-momentum transfer, the region of the zero sound branches can be crossed for 
some values of the velocity in the medium: If (uj, q) coincides with a zero-sound solution 
k% s = (io zs ,q zs ) of the dispersion relation D a _ LU L(k^ s ) = 0, there will be a critical value of the 
relative velocity v solution of 

lo zs = uj, e\ = jvpli: sin a cos 99 sin 6 — cos a(l =F cos 6)] (14) 

(c/. Eqs. (^,13)) for which the real part of the dispersion relation vanishes. Since the disper- 
sion relation enters in the denominator of the propagators (cf. Eq. (0)), a pole appears in 
the expression of the cross section (cf. Eqs. (§,^)) if the imaginary part of D^-uL vanishes 
as well. For velocities near v, D^^^l is still reduced and as a consequence, the cross section 
enhanced. 

For example, at n=riQ and T=0, we have drawn on Fig. 2 the path swept by the compo- 
nents of momentum transfer for fixed laboratory energy (7) a b=300 MeV) and relative velocity 
as the integration needed to calculate the total elastic in-medium cross section is carried out. 
It is seen that it falls on the upper part of the zero sound branch for v y =0Al corresponding 
to the location of the peak in Fig. 1. 
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The lower branch of the zero sound falls in the particle-hole damping zone. On the 
other hand, part of the upper branch falls outside the particle-hole damping zone where 
the imaginary part of D a _ LO i vanishes corresponding to undamped modes at our level of 
approximation. This is at the origin of the narrowest peaks shown in Fig. 1, while the 
smoother peak is to be related to the lower branch for which particle/hole damping is active. 
The zero-sound branches become more important as the density increases above a threshold 
value, their importance is reduced as temperature increases. It is indeed to be seen in Fig. la 
that the peak structure begins to develop above n=0.5no and in Fig. lb that it is washed 
out at high temperature. 



4 Discussion and conclusion 

In a relativistic OBE model of the N-N interaction with a, u and ir meson exchange, we 
investigated the observable consequences of a modification of the elastic proton-proton cross 
section inside symmetric nuclear matter by screening of the interaction, i.e. by approximating 
the meson propagators at the RPA level. In particular we concentrated on the influence of a 
non vanishing relative velocity between the center of mass of the collision and the rest frame 
of the fluid. The velocity dependence was found to be non monotonous. For moderate values 
of the velocity, the cross section is enhanced at low densities and reduced at high densities 
if the velocity is inside the reaction plane (y, z) and no appreciable modification is observed 
for velocities orthogonal to the reaction plane (v//x). 

The main result of our calculations is the enhancement introduced by the excitation of 
the zero sound modes of the medium on both the differential and the total elastic cross 
sections. A peak appears for a charateristic value of the relative velocity. Its location is 
density dependent; the peak is stronger at higher density and gradually blurs out when 
increasing temperature. It occurs for v//y at physically relevant values 0.3 < v y < 0.8. but 
not for v//x or for unrealistic values of vjjz. 

Such an asymmetric configuration has not yet been considered up to our knowledge. A 
selective enhancement of the cross section for a well defined value of the density and velocity 
will probably have a characteristic signature in terms of the transverse flow and angular 
distributions. In a naive hydrodynamical picture, we found that the collisions are more 
efficient when they take place in a background matter flowing in the direction corresponding to 
the "sidewards flow" , moreover the scattered products are predominantly directed backward 
and forward. However, it is delicate at this stage of our investigations to put forward a definite 
prediction on the consequences on experimental observables, since they are the result of the 
interplay of many factors. For example, it is known that viscous effects Q or the softening of 
the equation of state (EOS) through a possible phase transition to the quark gluon plasma 



[16 1 may reduce the directed sidewards flow. We think that our findings may bring new 
arguments in the context of the anisotropic flow observed recently both in Au+Au collisions 
at 100-800 MeV/c at the FOPI facility at GSI Darmstadt and by the E877 collaboration at 
AGS in Au+Au collisions at 10.8 GeV/c. The anisotropy is not yet fully understood fl7|l 
and has been analysed mainly in terms of the EOS. It could be worthwhile to investigate the 
effect of an anisotropically enhanced cross section in this kind of numerical simulations. 

The enhancement we have found only depends on the reduction of the real part of the 
dispersion relation associated with the vicinity to a zero sound mode in the matter. The 
existence of a zero-sound branch in nuclear matter is well documented (see e.g. [[H], |l^]), so 
that our result is a quite general one. Its size depends on the parameters, in particular on 
the cutoffs A CT and A^. Our parameters come from a fit of the free cross section; the cutoffs 
A ~ 1200 MeV and coupling constants are rather lower those of the Bonn potential, so that 
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we feel we are on the safe side. Let us say a few words on higher order corrections. Vertex 
functions such as calculated by Allendes and Serot [18| are amenable at low and intermediate 



momentum transfer to a monopole form factor with A ~ 1100 MeV. If the coupling constants 
are decreasing functions of the density as argued from Dirac-Brueckner-Hartree-Fock calcu- 
lations, the density dependence of our enhancement can be modified. (The determination of 
a more accurate density dependence would require a fully consistent calculation of ladders as 
well as loops which is not yet available.) Even if these branches would disappear completely 
(by applying very low cutoffs), there would still be a strong reduction of the dispersion re- 
lation D a - u L(u>,q) in this zone due to a — uj mixing and consequently an enhancement in 
the cross section. We believe that this broad enhancement at finite relative velocity matters 
more than the presence of a high but narrow peak. Its location will still be determined by 
the upper limit of the particle hole damping zone. 

Finally, we would like to warn the reader against a possible confusion with an enhancement 
found by Aim el al. Q in an other context. These authors perform a nonrelativistic Brueckner 
calculation. Their enhancement is caused by a pole in the propagator of the intermediate 
two-nucleon state; the result is interpreted as a precursor of superfluidity. The enhancement 
occurs for low values of the density ti/uq ~ 0.5; the peak is reduced and finally disappears 
when the total momentum K of the pair deviates from zero. This is clearly different from 
our case where the enhancement is due to a pole in the mixed meson propagator and occurs 
for high values of the density equal or larger than n/no = 0.5 and for nonvanishing values of 
the velocity v/c= 0.3 - 0.7. 

References 



[8] 
[9] 

[io; 
[ii 

[12 



D. T. Khoa et al, Nucl. Phys. A548 (1992) 102. 

E. I. Tanaka, et al, Phys. Rev. C52 (1995) 316. 

G.D. Westfall et al, Phys. Rev. Lett. 71 (1993) 1986; V. de la Mota et al, Phys. Rev. 
C46 (1992) 677; B.A. Li, Phys. Rev. C48 (1993) 2415. 

D. Klakow, G. Welke and W. Bauer, Phys. Rev. C48 (1993) 1982. 

see e.g. R.B. Clare and D. Strottman, Phys. Rep. 141 (1986) 177; H.H.K. Tang and 
C.Y. Wong, Phys. Rev. C21 (1980) 1846, ibid C26 (1982) 284; W. Schmidt et al, Phys. 
Rev. C47 (1993) 2782. 

A. Bohnet et al, Nucl. Phys. A494 (1989) 349; A. Faessler, Nucl. Phys. A495 (1989) 
103c. 

C. Fuchs, L. Sehn and H.H. Wolter, Nucl. Phys. A601 (1996) 505. 

G.Q. Li and R. Machleidt, Phys. Rev. C48 (1993) 1702. 

T. Aim, et al, Nucl. Phys. A587 (1995) 815 and Phys. Rev. C50 (1994) 31. 

B. ter Haar and R. Malfliet, Phys. Rev. C36 (1987) 1611. 
J. Diaz-Alonso, L. Mornas, Nucl. Phys. A629 (1998) 679. 

Nijmegen NN-Online Database, V.G.J Stoks et al. Phys. Rev. C47, (1993), 761, ibid. 
C48 (1993) 792 and SAID program by R.A. Arndt et al. Phys. Rev. C50, 2731 (1994). 



7 



[13] J. Diaz Alonso and A. Perez Canyellas, Nucl. Phys. A526 (1991) 623. 

[14] E. Gallego, J. Diaz Alonso and A. Perez, Nucl. Phys. A578 (1994) 542; J. Diaz Alonso, 
E. Gallego and A. Perez, Phys. Rev. Lett. 73 (1994) 2536. 

[15] J. Diaz Alonso, B. Friman, P.A. Henning, "Normal Modes in Nuclear Matter at Finite 
Temperature". Unpublished preprint G.S.I. Darmstadt (1992). Some results of this last 
reference have been published in P.A. Henning, Phys. Rep. 253 (1995), 235. 

[16] D.H. Rischke, Nucl. Phys. A610 (1996) 88c. 

[17] P. Crochet et al, FOPI collaboration, Nucl. Phys. A627 (1997) 522; J. Barrette et al, 
E877 collaboration, Phys.Rev. C55 (1997) 1420. 

[18] M.P Allendes and B.D. Serot, Phys. Rev. C45 (1992) 297. 

[19] A.L. Fetter and J.D. Walecka, "Quantum theory of many particle systems", McGraw- 
Hill, New York (1971) p. 183#.; J.C. Caillon, J. Labarsouque, Phys. Lett. B352 (1995) 
193. 

Figure captions 

Fig. 1: Ratio of the in-medium cross section with respect to its value in the vacuum 
&med/o~free as a function of velocity (0,^,0). — (a) for density n/no = 0.1, 0.5, 1, 2, 5 — (6) 
temperature dependence at n = uq. 

Fig. 2: Levels of constant value of the mixed a— longitudinal u> dispersion relation -D CT -wL 
in the q - u) plane of transferred momentum. The long-dashed line v y = 0.41 corresponds to 
the resonant condition *f y v y psm9 = uj zs . 
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Fig. 1a 
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Fig. 2 



contour plot of the mixed S-VL dispersion relation, T=0, n/n0=1 




